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Abstract. We present a calculation of the weak localisation correction to the 
scattering and to the transport mean free path, for waves propagating in a 5- 
correlated random potential, going beyond the usual diffusion approximation 
for the loops in the interfering path amplitudes. We find that not only the 
transport mean free path, but also the scattering mean free path is affected by 
the interference contributions. We also find the dependence of the leading order 
contribution to both, the scattering and the transport mean free path, on the 
disorder parameter l/k£ to be linear. This is in contrast to the result obtained 
from the diffusion approximation according to which the scattering mean free path 
remains unaffected when changing the disorder parameter, whereas the correction 
of the transport mean free path scales like l/{k£)^. 
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1. Introduction 

A thorough understanding of the propagation of waves in disordered potentials is of 
great importance in many areas of physics. Classically, one expects diffusive transport 
of the average wave intensity through the medium, i.e. one expects that interference 
effects play no role, on average. In other words, the neglect of interference results 
in the average wave intensity performing a random walk just like a classical particle 
would do. This picture was applied quite successfully to the propagation of electrons 
in conductors leading to the ohmic dependence of the transmission on the sample 
thickness [T] and, in optics, to the theory of radiative transfer [2]. 

In general, diffusive transport is characterized by the fact that the square of 
the distance r of the particle from its initial position increases, on average, linearly 
with time, i.e., (r^) ~ 2dDt, with diffusion constant D, in d dimensions. We will here 
consider the case c? = 3. If particles are continuously emitted from a source, the particle 
density evolves into a stationary distribution, which, in the diffusive case, decays 
asymptotically like l/r, with a prefactor which is solely determined by the strength of 
the source and the transport mean free path £t (see section 12.41 for more details) . The 
transport mean free path thus characterizes the behaviour of the stationary particle 
density at large distances r from the source, or - vice versa - the latter can be used to 
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define the former. On the other hand, the scattering mean free path tg, i.e., the average 
distance between two successive scattering events, determines the density of particles 
at small distances, where the particles have not yet been scattered after emission from 
the source (for waves, this particle density corresponds to the coherent intensity). 
This contribution to the density of particles decays exponentially on the length scale 
is- It is thus possible to measure the scattering mean free path via detection of the 
coherent intensity as a function of the sample thickness. For isotropic scattering of 
classical particles, transport and scattering mean free path are identical, whereas, in 
the general case, they are related by £t = 4/(1 — (cos(0))), where (cos(0)) is the 
angular average over the cosine of the scattering angle. The transport mean free path 
therefore defines the length scale over which the particles' direction of propagation 
is randomised. For charged particles such as electrons, another important quantity 
describing stationary transport is the electrical conductivity cr, defined by the diffusive 
current in response to an externally applied electric field. In the classical case, the 
conductivity is proportional to the transport mean free path, i.e., a = ne^lt/'mv, 
called Drude conductivity, with electronic density n, charge e, mass m and transport 
velocity v (see, e.g., chapter 7.2 of [T). Finally, also the diffusion constant is related 
to £t via D = Itv/i. 

Coming back to the transport of waves, the above considerations referring to the 
random walk of a classical particle neglect one important feature inherent to waves - 
their ability for interference. The Boltzmann approximation - which yields the above 
classical results, and which is valid for wave lengths much smaller than the scattering 
mean free path - assumes that, upon averaging over the disorder, all interference 
contributions vanish. However, if the system exhibits reciprocity symmetry [3] certain 
interference contributions do survive the disorder average and lead to deviations 
from the classical results. The most important contribution originates from the so 
called maximally crossed scattering paths |4j|5], where the two amplitudes ip and ip* 
forming the wave intensity visit the same scatterers, but in opposite order. In 
a backscattering experiment, these amplitudes will interfere constructively in retro 
reflection giving rise to the famous coherent backscattering effect [5]. Analogously, 
this effect enhances the probability for the wave energy to return to any point inside 
the medium and thus, because of flux conservation, reduces the overall transport over 
large distances - an effect known as weak localisation [7]. Nevertheless - in three 
dimensions and if the disorder is not too strong (i.e. if the scattering mean free 
path is still larger than the wavelength) - the propagation of the average intensity 
still resembles a diffusion process, although the two parameters is and it - which, as 
discussed above, govern the behaviour of the stationary intensity distribution on small 
and large length scales, respectively - are then modiflcd with respect to their classical 
values. In particular, the weak localisation effect manifests itself as a reduction of the 
transport mean free path. 

In most of the standard literature, e.g. [DIHIIS], the weak localisation correction 
is calculated for the diffusion constant instead of the transport mean free path, and, 
furthermore, a cut-off procedure is usually employed which relies on the following 
argument: as described above, weak localisation is related to the probability for 
the wave intensity to return to a position where it has been previously. To obtain 
a finite quantity, one considers a finite small volume around the initial position to 
determine the return probability. In momentum space, this amounts to a cut-off at 
large momenta the precise value of which depends on the model under consideration 
[lOl [Tl] . Furthermore, the return probability is usually evaluated within the diffusion 
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approximation, assuming that many scattering events occur before the wave intensity 
returns to its previous position. As a result of these approximations, the reduction 
of the diffusion constant turns out to scale like l/(fc£)^, where k denotes the wave 
number. 

In the present paper, we relax these approximations and calculate the corrections 
of the scattering and of the transport mean free path without employing any 
assumption apart from ki ^ 1 ('weak disorder'). For simplicity, we will concentrate on 
the case of a white-noise disorder potential, where, in the Boltzmann approximation, 
scattering is isotropic, and, hence, the classical values for the scattering and the 

(cl) fell 

transport mean free path are identical: il = Q — £■ We will then determine 
the leading modifications of is and it due to interference effects for ki ^ 1. For 
both quantities, the leading corrections turn out to scale like l/{k£). Let us note 
that a similar calculation of the leading weak disorder corrections has already been 
performed for the conductivity a (which is defined by the response to an external field, 
e.g. in the case of electronic transport) [121 US]- In the present paper, we repeat this 
analysis for the scattering and transport mean free path, which - as explained above 
- are general quantities determining the stationary intensity distribution for any kind 
of wave undergoing a diffusive scattering process in the absence of an external field 
(which is the typical situation, e.g. for light waves or sound waves). As we will see, 
the weak disorder correction of a is identical to the correction of the product of the 
scattering and of the transport mean free path. 

The paper is structured as follows: 

• The next section presents a short review of the diagrammatic calculation of the 
average wave intensity. For this, we first introduce the model under investigation 
and derive the vacuum Green's function. Then, the disorder average is introduced, 
to obtain the average Green's function. This enables us to define the scattering 
mean free path ig. With the help of the average Green's function we write down 
a Bethe-Salpeter equation for the average intensity propagator which is solved 
within the ladder approximation and describes incoherent (or classical) transport. 
We close the section with the derivation of a steady state diffusion equation for 
the average intensity on large length scales. The transport mean free path it is 
defined with the help of this equation. 

• In secti!07i[^interference corrections of the mean intensity are considered. We first 
specify the types of diagrams we include into our calculations, and show how these 
affect the transport mean free path it- We show that these additional diagrams 
also require a change of the scattering mean free path is. The evaluation of the 
corrections of is and it is sketched, and analytical results for the leading order 
terms are derived. 

• The results obtained in section |3] are summarised and complemented in section^ 
The leading order contributions to the corrections of the scattering and transport 
mean free path are discussed, and the results are related to previous studies on 
disordered electronic systems concerning the weak localisation corrections of the 
conductivity. 

• The last section concludes this article, and gives a short outlook on future 
perspectives. 
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2. Green's functions, disorder average and classical contribution 

After introducing the random wave equation as the starting point of our 
considerations, and giving the reader a short reminder about diagrammatic techniques 
for performing the disorder average, we will focus in this chapter on the regime of 
incoherent transport. Within the ladder approximation, valid for weak disorder, the 
average intensity is found to follow a classical random walk, with the average step 
length given by the scattering mean free path ig. On large length scales, the random 
walk reduces to a diffusion equation, which, as described in the introduction, defines 
the transport mean free path. 



2.1. Model 

We investigate the propagation of classical waves in a three dimensional disordered 
medium described by the stationary Helmholtz equation 

(l^ + e{l + V(r))y(r)^3{r), (1) 

where V' denotes the field, k the wave number, and j is the source term. The potential 
V , describing the disordered medium, is a random function which is assumed to be 
uncorrelated (white noise) and to exhibit Gaussian statistics: 

{V{r)) ^ 

47r (2) 
k^{V{r)V{r')) = —5{r-r'). 

Here, 47r/£ is a measure for the strength of the potential parametrised by a length 
scale In the limit of a very weakly disordered medium, the scattering mean free 
path turns out to be given by see section [2?3l below. We use the Green's function 
method to treat the problem. For the Helmholtz equation with V = i.e. the vacuum 
case, the retarded momentum space Green's function Go is well known to be of the 
form (see, e.g., [M], chapter XH.9) 

Gob) = ,2 V . , (3) 

with 77 an infinitesimal positive quantity. In real space, the retarded Green's function 
then takes the form of a spherical wave 

Go(0 = , (4) 

47rr 

and the advanced Green's function is the complex conjugate Gq. 

In the presence of a scattering potential V , the Green's function g can be expressed 
in terms of the vacuum Green's function Gq and the potential V by means of the well 
known Born series 



g(r,0 =Go(r,r') + j driGo(r, ri)^(ri)Go(ri, r') 



+ j dridr2Go(r,ri)T/(ri)Go(ri,r2)F(r2)Go(r2,r') + ---, 
which can be recast into a Lippmann-Schwinger equation: 

g{r, r') = Go(r, r') + j driGo(r, ri)V^(ri).9(ri, r'). (6) 



(5) 
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To shorten the notation one expresses integral equations of this type with the help of 
diagrams: different types of lines represent the different types of Green's functions 
{Go{r), g{r,r') or, later, the average Green's function G(r)), crosses denote the 
scattering potential. In the next section we will perform the averaging over the disorder 
potential V. The appearing correlators {V{r)V{r')) are represented diagrammatically 
by crosses connected via a dotted line. As opposed to the Green's function g for 
a specific realisation of the disorder potential, the average Green's function G is 
translationally invariant. 



2.2. Disorder average 

The average Green's function G{\r — r'\) in presence of the random medium can be 
expressed via a Dyson equation (see, e.g., |14| . chapter XIII. 23): 

G(|r - r'\) = Go(|r - r'|)+ J dr,dr^Go{\r - r,m\r, - r^DGUr^ ~ r'|), (7) 

where S(r) is the self energy defined as the sum of all irreducible diagrams, i.e. 
diagrams which do not fall apart upon cutting a single Green's function, see figure [T] 
From (O, we obtain an algebraic solution for the Fourier transform of the average 
Green's function G, in terms of the vacuum Green's function Go and of the self energy 
S: 

G{p) = ^ ^ ^ ^ = ^—^ — . (8) 

G^\p)-j:{p) k^-p^^^p) 

Here, we do not need the infinitesimal displacement of the poles rj like in ([3]), since 
the self-energy already possesses a finite imaginary part (see below) which shifts the 
poles away from the real axis. For weak disorder, i.e. fc£ 3> 1, the self energy depends 
only weakly on the momentum p and we can therefore write: 

G{p) = , (9) 

fc^ — p'^ 



with the effective wave-vector k — ^Jk'^ ~ I](fc). Neglecting the real part of the self 
energy (which can be absorbed by a suitable redefinition of the wave number fc), this 
can be written as fc = fc + i/2£s, with the scattering mean free path 

4 = ^ . (10) 

ImS(fc) ^ ^ 

In real space, the retarded Green's function takes the form 

G[\r-A)^- ^''f\'-:'\ (11) 
47r |r — r'l 

and the advanced Green's function is again the complex conjugate of pT|) . We see that, 
due to the disorder, the real space average Green's function decays exponentially on a 
length scale given by twice the scattering mean free path fg. Therefore, we will later 
determine the scattering mean free path as the decay length of the squared absolute 
value of the average Green's function, see In the next section we will see that Is 
is the length scale on which the coherent intensity decays. 

To investigate the mean intensity I = (ipip*) (the brackets denote the average over 
disorder realisations), one defines the average intensity propagator ^ as the disorder 
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S = 



+ 



U 



Figure 1. Diagrammatic representation of the first terms of the self energy E 
and of the irreducible vertex U. Crosses connected by a dotted line represent the 
two-point correlation function {V{r)V{r')) of the random potential, see ((2]|. Solid 
lines correspond to vacuum Green's functions (in E) or average Green's functions 
(in U). 



averaged product of a retarded with an advanced Green's function. The average 
intensity propagator fulfils the Bethe-Salpeter equation [I] 

$(ri, r2; r,, r,) = G{\r, - r:,\)G*{\r^ - r,\) 

+ j dr^ ■ • •dr8$(ri,r2;r5,r6)J7(r5,r6;r7,rg)G(|r7 - r-3|)G*(|r8 - r^W 
where U denotes the irreducible vertex or intensity operator, see figure [TJ 



2.3. Weak disorder approximation 

In the case of weak disorder, i.e., ki ^ 1, one may apply the so-called ladder 
approximation. Here, only the first term of the self-energy and the irreducible vertex 
are considered (see figure [Ij. In this approximation, the self-energy and the irreducible 
vertex read 



^{\r-r'\)^^d{r~r')Go{\r~r'\) 
An 

U{r5,re;rj,rs) ~ tH^s ~ re)6{r5 ~ r7)(5(r6 - rg). 



(13) 
(14) 



Together with lfTO|) . this approximation yields = ^- Using U of the form ([T4|) in 
the Bethe-Salpeter equation (IT^ . and assuming a (5-source of strength sq at Tq, i.e. 
jX*") — \^S{r — ro) in ([T]), then leads to the following transport equation for the 
mean intensity I{r) = $(ro, Tq; r, r): 



I{r)=Io{r) + J dr'Po{\r-r'\)I{r'), 



(15) 



where Io('") = |G(|ro — r|)|^ denotes the coherent intensity, i.e., the intensity 
arriving at point r without being scattered, and 



Po(|r-r-|)^^|G(|r-.-|)|-^ -P("';-y) 
*• 47rl; 1 — r'\ 



(16) 



Although is equals £ within the ladder approximation, we still keep is in the 
exponential since, later, when additional interference diagrams are taken into account, 
this equality will not hold any longer and a distinction between is and i is crucial. As 
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mentioned in the introductory section [1] the coherent intensity decays exponentially 
on the scale of the scattering mean free path, which can be seen from PT|) together 
with PH)) . 

The transport equation (jl5p can be interpreted in terms of a classical random 
walk. According to dTSl) . this random walk is characterized by the source term Io{r) 
and the probability distribution Po(|r — r'|) of the step length |r — r'| between two 
successive scattering events at r and r' . Remember that this simple form of the 
transport equation only applies to isotropic scattering from a ^-correlated potential. 
In the case of anisotropic scattering (i.e., if the random potential exhibits correlations 
on a length scale comparable to the wavelength), one can derive, instead of (flS)) . 
a radiative transfer equation for the intensity I(r, s) at point r, with propagation 
direction s fW. 

The solution of the transport equation (jlSp can be expressed in terms of the 
ladder propagator ^d^"! — '"2I), which we will need as a building block for the weak 
localisation diagrams in section [3] 

I(r)= y"dridr2|G(|r-r2|)|'L(|ri-r2|)Xo(ri). (17) 

Equivalently to ()15|) . the transport equation for the ladder propagator reads: 

47r f 
L{\ri - rsl) = —S{n - ra) + J drgLdn - rgDPodra - rsl) , (18) 

see figure [5J In Fourier space, the sequence of single steps reduces to a geometric 
series: 

Lil)--T.{Poil)) a.ctan(,4 T ' ^^^^ 

n=0 ^ e q£^ 

where we used the Fourier transform Pq of ([T6l) 



~ = 4 arctan(^ 

£ q£s 

Summing up all ladder diagrams without the first (N — 1) terms leads to the ladder 
propagator with at least N scattering events: 

7 / ^ 47r ^ /"S/nA" 47r (f arctan(g4)M)^~^ 
Thus, Li = L holds. 

Similarly, we define the crossed propagator C as the sum of all maximally crossed 
diagrams [5] (see figure [2]), for later use. By reciprocity symmetry [3j, C equals 
^2- Furthermore, we can construct propagators Cn consisting of maximally crossed 
diagrams without the first (N — 1) terms. It then holds that Ci = C, and 

CN{q) = Ln+i{q) N = l,2,3,..., (22) 

generally. Diagrammatic representations of the propagators Pq, L, L2 and C are given 
in figure m 
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2.4- Diffusion approximation & transport mean free path 

If we consider transport on large length scales (compared to £s), we may derive a steady 
state diffusion equation starting from (jl5l) . For this, we perform a Taylor expansion 
of the mean intensity X(r') around position r inside the integral in ([15]): 



I{r + r') = I{r) + r' • VZ(r) + K^r^ + . . • , (23) 

where summation over repeated indices is implied. Stopping the expansion at the 
second order and insertion into ()15|) yields 

I{r) -Io(r)+ j dr'Po(r') (l{r) + r' ■ VZ(r) + Ky^^I^^ .(24) 

The first term in the integral cancels the left hand side of the equation, because of the 
normalisation of Pq : 

J drPoir) = 1 . (25) 

The second term vanishes because it is an odd function of the integration coordinates. 
Thus, only the term quadratic in r' survives inside the integral: 

0=I,ir) + lJdr'P,ir')ry^^^. (26) 
By exploiting the spherical symmetry of Pq, we can use the identity 

y dr'Po(r')rf = y dr'Po(r>f , Vz,j (27) 
and conclude that 

J dr'Poir') ry^^^ ^\ j dr' P,{r'y^ M{r). (28) 

where A = denotes the Laplace operator. Thus ([^5)) takes the form of a stationary 
diffusion equation: 

I?(")AI(r) +Io(r) = 0. (29) 

It only depends on a single parameter 2?^^^ which is the stationary analogue of 
the diffusion coefficient appearing in the time dependent diffusion equation. The 
parameter is given by the relation 

= drr2Po(r-). (30) 

By comparison with the fundamental steady state diffusion equation, found in jl5J 
chapter 9, we see that 

^ ^ (31) 

3 

which we use to calculate the transport mean free path. The interpretation of Pq 
as the step length distribution (cf. section 12.31 above) then means that we calculate 
the transport mean free path via the variance of the step length distribution or, 
equivalently, by the mean square displacement. 

Within the ladder approximation, we find from (|25l). together with ([T6|) . 



1 ^'^ [ A \ni w"^ [a exp(-r/4) 4 

l = -jdr\G{r)\ =ydr^-^ = -, (32) 



Weak disorder corrections of the scattering and transport mean free path 



C 







X 

i H 


h n - 


h TTT + 






X 


x-x 






X 



LA = 



x-x 



x..><, x.x 
x'")< V'x 









L2 





Figure 2. Diagrammatic representation of the "single step" Pq, the ladder and 
crossed propagators L, L2 and C, and of the Bethe-Salpeter equation I II8I I satisfied 
by the ladder propagator. By reciprocity symmetry, C and L2 have the same 
value, see l|22| l. 



and thus = ^, which is consistent with the result found from using the definition 
(fTO]). Furthermore, inserting (HH) into (l30|) yields P^*') = £^/3 = £'^/^, Combining this 
with (PT|) we obtain: 

it=i = is. (33) 

Again, remember that this is only true for isotropic scattering - otherwise 
the scattering and transport mean free path would differ already in the ladder 
approximation . 

To close this section, it is instructive to look at the solution of (l29l) for a 6- 
source of strength sq, i.e., loif) = sq \G{\r\)\ = sq exp(— r/4)/(47rr)^ in an infinite 
medium. For r ^ ig, the exponentially localized coherent intensity can be replaced 
by Xo(r) = soigS (r) /47r, for which the solution of ([29|) reads 

3so 



I(r) 



(34) 



167r2£tr ' 

This is the asymptotic 1 /r decay mentioned in section [TJ with the proportionality 
constant 3so/167r^£t thus defining the transport mean free path. For given source 
strength sq, a measurement of the transport mean free path can therefore be performed 
via the detection of the average intensity at a large distance r from the source (where 
the coherent contribution is sufficiently attenuated). 



3. Weak disorder corrections 



In the previous section, we saw that, for isotropic scattering and within the ladder 
approximation, scattering and transport mean free path are equal, (I33|) . However, in 
the presence of reciprocity symmetry, additional diagrams, apart from the ladders, 
play a role. These contributions will give rise to corrections of the mean free paths 
and will therefore lead to deviations from (1551) . 

In this section, we will introduce these additional diagrams and show how they 
are incorporated into the transport equation for the average intensity. We will then 
calculate the corrections due to these diagrams, first for the scattering and then for 
the transport mean free path, and identify the leading contributions in the disorder 
parameter l/{k£). 
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Figure 3. Diagrammatic representations of the weak localisation contributions 
to the step length distribution P(|r — r'|). Each diagram is depicted in two 
different but equivalent representations: In the upper diagrams, the upper (lower) 
horizontal lines denote retarded (advanced) Green's functions, crosses connected 
by dashed lines refer to the two-point correlation function of the potential, see 
((2]l - which, in turn, amounts to a scattering event - and L2 and C represent the 
ladder and crossed propagators defined in ((2]l. In the lower type diagrams, solid 
(or dashed) lines denote average Green's functions G (or G*), full circles refer to 
scattering events, open circles to a specific point in space (without scattering 
event), whereas paired solid/dashed lines, (■■■■■•) or (•™"), indicate ladder or 
crossed propagators, respectively. 



3.1. Contributing diagrams 

We start by specifying the types of diagrams we include into the irreducible vertex 
U . Obviously, we have to consider the ladder contribution, see (|14p . which accounts 
for classical, i.e. interference- free, propagation. Interference-free means that the two 
amplitudes which constitute the intensity are perfectly in phase at any instance of 
the scattering sequence. Since, in the presence of reciprocity symmetry, the two 
amplitudes forming a maximally crossed diagram do not collect any phase difference 
either, these have to be included as well, see (P^ . When a maximally crossed diagram 
is connected with ladder diagrams on both sides, a crossing is generated. This crossing 
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introduces a phase difference between the involved amphtudes and thus vanishes in 
the limit of k£ — oo, i.e., the corresponding diagram will yield a small contribution 
in the weak disorder limit. A crossing can be decorated with one two-point potential 
correlation function {V{r)V{r')) i.e. crosses connected with a dotted line, without 
rendering the diagram higher order in the disorder parameter l/k£. In the case of the 
maximally crossed diagrams, these crossings are known in the literature as four point 
Hikami-boxes pLS,. Similar structures appear if one additional correlation function is 
introduced inside ladder diagrams, cf. figure |3]d) - i), whereas adding more than one 
of these correlation functions leads to contributions of higher order [TJ [T3] . Since any 
of these irreducible vertices, depicted in figure [3l may be inserted between two ladder 
scattering events at r and r', they will lead to a modified step length distribution 
P{\r — r'l) (see below). The diagrammatic representation in figure [3] is given in two 
different, but equivalent ways: the upper one is already familiar from figure [1] and 
figure [21 whereas the lower one corresponds to a real-space representation emphasizing 
especially the loop- like character of the additional diagrams. The diagrams carrying an 
index H are the familiar Hikami-boxes |16| . Since our aim in this article is to calculate 
the leading order weak localisation corrections, we do not consider a whole group of 
"small" diagrams (e.g. the diagrams given by the last two terms, the eighth and the 
ninth, of the irreducible vertex depicted in figure[T]) which contribute in order l/(fc£)^ 
and which can be found in [T3]. With this selection of the contributing diagrams 
we completely account for the corrections in l/k£ and ln(H)/(fc£)^ completely, see 
section |4l The mean intensity can be constructed from all possible combinations of 
these building blocks, together with the ladder building block Pq (figure [2]). This is 
achieved with a Bethe-Salpeter equation like (|15p 



where, now, in comparison to (I15|) . Pq is replaced by the modified step length 
distribution P which, apart from Pq, contains the additional diagrams shown in 
figure 131 see ([36]) below. 

Our approach is very similar to the one used in [17] , where the weak localisation 
corrections of spectral statistical properties, namely the two level correlation function, 
of two-dimensional disordered conductors were calculated for systems with spatial 
extent of the order of the scattering mean free path (denoted as "ballistic limit", 
there). These authors also use the exact propagator within the loops and the exact 
form of the Hikami boxes. In |17| . the model considered is a non-interacting electron 
gas on a two dimensional torus. The use of a finite medium is necessary, since in a 
two dimensional infinite medium all states are localised. Thus, the authors study the 
behaviour of the weak localisation corrections with respect to changes of the medium 
size. In contrast, we perform our calculations for a three dimensional infinite medium 
below the mobility edge, i.e. in the weak localisation regime. We analyse the weak 
localisation corrections of transport quantities, like the transport mean free path, 
under changes of the disorder strength. In the view of |17j . our approach is to allow 
the loops, which lead to the weak localisation corrections by influencing the return 
probability, to also probe short scattering paths - meaning the loops not only explore 
the "diffusive" , but also the "ballistic" scale of the medium. 

Lets return to the modified step length distribution P. Summing up all diagrams 
shown in figure [3] yields 




(35) 



P{r) ^ Po{r) + Pnir) + Peir), 



(36) 
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where 

Ph = Ph^ + PHn + Phc (37) 

and 

Pb = Pba + Pbb + Pbc + Pbd + Pbe + Pbf ■ (38) 

The respective diagrams fulfil the following symmetry relations: mirroring vertically 
is equivalent to complex conjugation (interchange of the upper line, G, and the lower 
line, G*), whereas mirroring horizontally (either the complete diagram or just the 
upper or lower half) does not change the value of the diagram (reciprocity symmetry). 
Thereby, we can infer 

Phb ^ P;i,, = Pb, = P*Bb, PBa^P*B^, Pbe=P*b,, (39) 

i.e. only four different diagrams have to be calculated: Ph ai Ph^i Pbc ^^id Pbe- 

As already pointed out before, the inclusion of additional diagrams into the step 
length distribution P entails additional contributions to the irreducible vertex U as 
compared to the ladder approximation. Since the irreducible vertex and the self- 
energy are linked via a Ward identity which ensures conservation of the total energy 
flux [T5], this also implies the inclusion of additional diagrams, apart from (fT5)) . into 
the self-energy. This results in a change of the scattering mean free path. In the case 
of (j35p . energy flux conservation is ensured by the normalisation of the modified step 
length distribution 

j drP{r) = 1, (40) 

from which we will calculate £s as shown below. As we have checked, the same result 
is obtained when calculating is directly from the self-energy, see ([TU]), using a Ward 
identity [TH]. On the other hand, the variance of the step length distribution yields 
the transport mean free path via 

i I dr r'Pir) = ^ , (41) 

as discussed above, cf. dSj). 

Finally, we give an example of how to translate the diagrams depicted in figure [3] 
into formulas. The diagram Ob), e.g., represents the equation 

PHA\r~r'\)^(^^^ Jdr,dr,dr,Gi\r~rs\)G{\r3-r,\)G*i\r~r2\) ^^^^ 

Ci\r, - r2|)G(|r2 - r,\)Gi\r3 - r'\)G*{\r, - r'\). 

We recognize the four Green's functions G, the two complex conjugate Green's 
functions G* , and the crossed propagator C. According to its definition, see (|22|) 
and figure [U C already contains two scattering events at the starting and end points 
(ri and r2, respectively). The remaining two scattering events (at and r') are 
taken into account by the prefactor (47r/£)^ originating from the potential correlation 
function We must not include a factor of Air/i for a scattering event at r, since the 
transport equation (I35p would otherwise double count these events upon concatenating 
the building blocks. 
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3.2. The scattering mean free path 

Since, according to pT|) . gives the decay length of the average Green's function, it 
follows from (ITCl) that 



J drPoir) - I 



(43) 



The scattering mean free path is thus obtained from the normahsation (^0)) of the 
total step length distribution P = Pq + Ph + Pb, (p6| . as follows: 



1 = 1-/ dr (PH(r) + Psir)) ^ 1 + ^ + ^ , (44) 

where we introduced the corrections of the scattering mean free path due to the 
different Hikami, SHg, and dressed ladder boxes, SBg. To evaluate the integrals, we 
switch to Fourier space, since there we know the analytic form of the ladder /crossed 
propagator L/C, see (PT|) and (P^. The corrections thus take the form (beware the 
minus sign from (|44l) !) 

= - / drdr'dri2 H{r, ri, r', r2)C(ri2) 

= - J dri2H{ri2)C{ri2) (45) 
:Hiq)Ciq). 



(27r)3 

Here, H{r,ri,r',r2) represents the sum of diagrams a)-c) (figure [3|) without the 
crossed propagators connecting ri and r2. As an example, the contribution from b) 
reads: 

JfB(r,ri,r',r2)- (^y) J dr,G{\r - r^min - n\)G* {\r ^ r,\) ^^^^ 

G(|r2-r3|)G(|r3-r'|)G*(|n-r'|). 

Integration over the outer points r and r' yields the reduced Hikami-box H{ri, = 
H(r 12) with Fourier transform H{q) = Hp^{q) + H-Q{q) + Hc{q). Inserting this into 
yields, correspondingly, SHg = SH a,s + SHb.s + SHc.s- In the same way, we get 

for the dressed ladder boxes Ba d (the contributions from Be.f are special and we 

consider them separately in a moment): 

SBa+...+u,s f dq 



{2n) 



B{q)L2iq), (47) 



with the Fourier transform of the dressed ladder boxes B{q) — BA{q) + • • • + B£,{q). 
The Fourier transforms of the boxes are readily obtained to be 

HA(q) = — ^( 2arctan(g4)-arctan((j4+2fc4)-arctan((74-2fc4) (48a) 



2 



i?B(g) = [^^j (^-2iarctan(g£) 

- ln(l - \{ql - 2ke)) + ln(l + i(g£ + 2H))) ' 
= H*ciq)=BAiq)=B*B{q), 



(48 &) 
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and, finally, 



Bciq) 



1 

V4fcg^ 

X [2i arctan((7^) + ln(-i + 2kl - qi) - ln(-i + 2kt + qi)] (48c 
X [-2i arctan((7£) + ln(i + 2kl - ql) - ln(i + 2kl + ql)] 

^Bu{q). 

Here, we keep the distinction between is and £ only in the _ffA-box (j48ap . whereas we 
set is = £ in the remaining boxes i?B,c and i?A,B,c,D- As we will see later, the latter 
boxes lead to higher order corrections as compared to Ha, such that the difference 
between is and i can be neglected. Furthermore, we used the symmetry properties 
discussed earlier, leading to Hb = Hq — Ba — B^, see ([39|) . The final integrations 
over q, including the propagators L2 and C, see and (H71) . have to be performed 
numerically. The results will be presented in sectional 

The corrections due to the diagrams from figure 12] h) and i) which do not 
contain full ladder or crossed propagators can be calculated directly (without going to 
momentum space), and analytically. Because of translational symmetry, we set r = 
and obtain (ri/2 = |'"i/2|) 



SBv. s /47r^ ^ 



£ 



J dr'dndra G(ri)G(|ri - r2|)G*(r2) 



G(|r2-ri|)G(|ri-r'|)G*(|r2-r'|) 



, ki \ , , -l + 2iki 
2k£ arctanh — + (i + ki) In 



2(H)2 
fSBp 



i + k£j \{i + kef 



(49) 



For ki 1, the sum of these two contributions has the asymptotic behaviour 
SBe+f,s _ SBe,s + SBf^s -1 + ln(2) - \n{ki) 



(50) 



i i /c£»l (H)2 

Moreover, we can obtain the leading correction of the scattering mean free path 
analytically. It turns out to be given by a part of diagram [3] a) , depicted in figure [H 
where only the first term {AnPo/i) in the scries forming the crossed propagator is 
taken into account: 

i^ aictanikis) 41n(l + (fc4)2) tt In (H) 

~ kP ^ 2fc2£3 ki^i ~2M ^ [kif ■ 

Again, we set is = i io obtain the two leading contributions on the right hand side 
of (jSTj) . since £s — > ^ for fc£ — 00. (The distinction between is and i will affect only 
the next leading contribution proportional to l/{kiY, see the coefficient 02 of 5Ha,i 
in table [U below.) 
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Figure 4. Short-loop diagram responsible for the leading correction of the 
scattering and transport mean free paths, depicted in two different but equivalent 
representations (as in figure|3]l. It is obtained as a part of diagram [3] a), where the 
loop between ri and r2 is replaced by the direct connection, or, in other words, 
only the first term in the series for the crossed propagator C, see figure [2] is taken 
into account. 



The remaining contributions to diagram |3] a) are given by replacing AttPq / £ with 
C2 — L3, see (1^^ . in psj) . i.e. by subtracting a single scattering event from the full 
C: 

^"^^'"^ - ^ M<l)C,{q) = I -P-.HaA<DC{<D (52) 



£ J (27r)3 J (27r)3 

where we used C2{q) — PQ{q)C{q), see (f2T|) and (|22l) . and introduced the extended 
ffA-box 

HA,2iQ)^Po{q)HA{q). (53) 
3. 3. The transport mean free path 

To calculate the corrections of the transport mean free path, we use (01]), which we 
rearrange in the form 

7^2^/ d^'^'(^o(r) + PH(r) + PB(r)). (54) 
Using (fT6|) we can evaluate the integral over the ladder step Pq which yields 



£ \ £ 2£J 



J drr^(^PH{r)+PBir) 



(55) 

£ J £ £ ' 

where, analogously to the case of the scattering mean free path, we introduced the 
corrections due to the Hikami- and dressed ladder diagrams, SHt and 6Bt, respectively. 
We can put (|55p into a more convenient form by writing £s = £ + d£s, and expand the 
square of £s/£ for small 6£s to obtain 

£t _ Js 5Ht 5Bt 

5H^ 6Bs SHt SBt ^ ' 

^ 1 + 2 — ^ + 2 — ^ H H -. 

£ £ £ £ 

We see that the corrections SHs and SBg of the scattering mean free path calculated 
in the previous section also enter into the corrections of the transport mean free path. 

In order to use the exact ladder/crossed propagator in the calculation of the 
corrections, we again switch to the Fourier transforms and obtain: 

f -2^/(l^^^^'(^)^(^) 
SBA+-+u,t 1 f dq (2)/ ^T ( N 
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where we introduced the Fourier transforms H^'^\q) and B^'^^{q) of the functions: 

H^'^\ri2)^ J drdr' \r -r'f H{r,ri,r',r2) (58) 



and 



B^^\ri2)= J drdr' \r -r'f B{r,ri,r',r2), (59) 

respectively. As compared to H{ri2) and B{ri2) in section [321 the additional term 
|r — r'l^ appears in ([SS)) and (|59|) . This makes the calculation more complicated, 



however - with some effort - the integrals in ([55)1 and ([5^ can still be performed 
analytically. 

(2) 

We start with the contribution H)^ . Using translational invariance, we choose 
the point r2 in diagram [3] a) as the origin. We then introduce the Fourier transforms 
of the Green's functions 



(2,) 

exp[ipi ■ r — ip2 ■ (r — ri) ^ ip^ ■ r' — ip4 • (ri — r')] ^^^^ 

G{pi)G*{p2)G{p3)G*{pi). 

Substituting the variables u = {r + r')/2 and v = r — r', shifting p2 by —P3, and 
using the relation 

J dr exp(ip • r)/(p) = ^{27r)H{p)Apf{p), (61) 

valid for sufficiently smooth functions /, we arrive at the expression 

„(2). ^ 47r f dp3dp4 . 
HX'{ri)=-— J -^^-^exp[-i(p3-p4) TiJ 

G{p3)G*{p4) Ap, [g{P2+P4)G*{P2+P3) 

The Laplacian of the product of Green's functions follows from ([S]): 
Ap, (G{p2+Pi)G*{p2+P3)) - 6&{Pi)G*{p3) + 6G(p4)G*2(p4) 

V / P2=0 

+ 8plG^p,)G*{p3) + 8plG{p,)G*^p3) (63) 
+ 8p4-P3G^(P4)G*2(p3). 

Choosing spherical coordinates, one can first perform the integration over the angular 
variables, and then evaluate the remaining integrals over p3 and pi using the calculus 
of residues. Finally, Fourier transforming the result with respect to ri yields the 
desired H^^ {q). From this, the correction is obtained by numerical integration of 
H'^\q)C{q) over q (cf. §7^). 

Like in section [321 the short loop (diagram in figure [4]) accounts for the leading 
order correction and can be obtained analytically again by replacing the full crossed 



(62) 

P2=0 ' 
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propagator C with a single step iirPo/i: 

SHA,i,t 1 /■ _dq ~ (2) 47r ~ 



1 



^''^f+\ktr^' + (^6^^^ + arctan(fc4) (54) 



+ (5-12(fc4)2)ln(l + (fc4)2 



5tt ln(M) 



The remaining contribution to SH^.t, i-e. 5-ffA,2.t = 6HA,t — SHA.i,t^ reads: 
again with the extended _ff^^2"box given by, of. (|53|) . 

i?S('z) = i^o(9)i?f (g). (66) 

Like SH^ i t, the contribution (Ji^E+F.t from diagrams [3] h) and i) can also be 
calculated analytically: 

5i?E+F,t 1 



2{k£y + 52{key + 32{key 



e ~ 24(fc^)4 [1 + 5(M)2 + 4(H)4J 

+4H [27 + 115(H)2 + %{Uf ~ m{klf] arctan(fc^) 
-20fc£ [3 + ll(kif - 8ik£)^ - 16(H)^] arctan(2fc£) 
- (10 - 46(H)2 - 440(fc^)'* - 384{kif) ln(l + (H)^) (67) 

+ (5 - 35(H)2 - 280(H)4 - 2A0{k£f) ln(l + 4(M)2 

31n(M) 2-51n(2) 
feC» 1 (H)2 ^ (fc^^p ■ 

For the contributions from the remaining diagrams, [3]b)-g), one can proceed 

analogously to the calculations Since the resulting expressions are very 

lengthy, we do not give them explicitely but rather discuss their general behaviour. 

— (2) — (2) — (2) — (2) 

A plot of the Fourier transforms of the functions Hj^ ; ^b+c ~ ^a+b ^^'^ ^c+d ' 

functions of momentum q, is given in figure[5j These functions are directly comparable 

with each other, since their contributions to the correction of the transport mean free 

path are given by integrating these functions multiplied by C{q) — ^2(9) over the 

momentum q, see ([57)) and (|65p . The inset of figure [S] shows that all functions remain 

finite for vanishing q. For these limiting values we obtain, in leading order of l/k£: 

i7g(0) « (^)' (68a) 

^4'ic(0)-4'|B(0)«-2(^)' (68&) 
B^cUo) « ■ (68c) 
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3 ■ 10" 



2 ■ 10" 



1 ■ 10" 



-1 ■ 10" 



-2 • 10" 



-3 ■ 10" 




Figure 5. The functions H'j^\{q) (solid), hI^Im 



250 



Sa|b(9) (dashed), 



and BQ^-^lq) (dotted), from the end of section 13.31 for the disorder parameter 
k£ = 100. Due to the very lengthy nature of the analytical expressions of these 
functions they are not given explicitely in the text. The inset shows that the 
functions remain finite for vanishing q. 



From figure [5l we also see that the function (g) is cut off at q 



2fc, whereas the 



functions -^b+c (?) ' ^a+b ^iid -Bc+d (?) peaked at ql = \/l + {2k£) 

and then decay like 1/q, asymptotically. 

In order to verify the consistency of our results, the values at g = given above 
can also be derived from the ordinary boxes -ffA,B,c and i3A,...,D, see (|48ap . (|485|) . and 
^cj), by using the following identity: 



J'{r'f{r)}{q) 



9=0 



-Aq/(g) 



J'il) 



9=0 



f'iq) „ , (69) 



valid for the Fourier transformation of a spherically symmetric, twice difFerentiable 
function / multiplied by r^. This relation can be derived with the help of (pT|) . 



4. Results 

In the previous section, we have derived expressions for the different contributions to 
the corrections of the scattering and the transport mean free path. Except for the cases 
of the special diagrams shown in figures HI [3] h) and i), which are given analytically, 
the integrals in these expressions have not yet been evaluated, cf. ([15|) . (|Tf)l . ([5^ . 
(|57p and ([55)1 . We now analyse the asymptotic behaviour of these contributions, by 
considering the data obtained from numerical integration. 



4.1. The scattering mean free path 

The leading asymptotic behaviour of (|5T|) . i.e. the correction due to the diagram 
depicted in figure IH is SHa,i,s/^ — > — 7r/2H for k£ oo. In higher order, it also 
contributes corrections of the order of In {k£) / (ki)'^, 1/ {ktf' , and further. 
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Table 1. Coefficients ai and 02 defining the corrections of the scattering 
mean free path in orders ln{k£) / {k£)^ and l/{k£)^, see 1701 1. for the different 
contributions from the diagrams in figure \3\ The data are obtained either 
analytically from l50t and 1511 or by fitting from list . l l47t and ((52J. 



SHA,l,s/i 




5HB+c.s/i 










= SBA+B,je 






ai 1 


-1.234 


1.234 


-1.234 


-1 


a2 l + 


-1.206 


0.155 


-2.258 


-1 + ln(2) 



The remaining diagrams give corrections at least of the order of hr(M) / (M)'^ 
and 1/ (M)^ or higher. We thus write: 

6H, SBs ai ln(H) + 

— —^ — m — f-fc^^-- (70) 

Table [T]hsts the coefficients ai and 02 in ([70)) for the individual contributions obtained 
analytically from ([50)) and (|5T|) (for (5i3E+F,s) or by fitting the numerically integrated 
contributions (cf. (|45|) and (|47l) '). We see that the different logarithmic contributions 
cancel each other (cf. the leading coefficient ai in table [1]), which only leaves 
contributions of the order of l/^ktf'. However, due to the contribution of SHa,i,s 
(|5ip. a term of the order of 1/ki, which thus constitutes the leading contribution, 
remains. The total correction of the scattering mean free path then reads 

T"^"2fc£+^^^^^^"'^ for 00. (71) 

This is the same correction one would find with the help of (jlOp from the second term 
of the self-energy depicted in figure [U see e.g. chapter 3.2 exercise 3.8 equation (3.82) 
in [T]. As already mentioned, this leading contribution originates from the short- loop 
diagram shown in figure HI where only a single step AttPq /£ is inserted instead of the 
full propagator C in the Hikami box Ha, see ([^T|) . As we have checked, inserting two 
steps 4:TtPq/£ into Ha - or, equivalently, a single step AttPq/I into the extended box 
Ha,2, see (|52l) - reproduces the next order correction, i.e., the coefficient ai — —1.234 
of 6Ha,2,s/^- Likewise, the ai-coefficients of SHb+c.s/£, SBa+b,s/£ and SBc+d.s/^- 
(which, as mentioned above, cancel each other together with SHA,2.s/i), are obtained 
by inserting a single step into the boxes Hb+c, Ba+b, and _Bc+d, respectively, i.e., 
(|45p and (|T7|) with AttPq/I instead of C = ^2- Therefore, all contributions from 
longer loops between ri and r2 (see figure [3]), i.e., those obtained by inserting a fuU 
propagator C2 = with at least three steps into the boxes Ha^2, Hb+c, Ba+b, and 
Bc+u, yield corrections of the order l/{k£)'^. As we have checked, the a2-coefficients of 
these longer loops also cancel each other. The same remains true if these longer loops 
are approximated by a diffusion propagator (proportional to 1/9^). For this reason, 
no correction of the scattering mean free path is found if the return probability is 
evaluated within the diffusion approximation, as mentioned in the introduction. In 
contrast, our exact calculation indicates a non vanishing correction of ig due to short 
loops. Let us note that, in a similar context, the infiuence of short scattering paths 
on weak localisation in quantum billiards has recently been analysed in |201 121j . 
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Table 2. Coefficients 61 and 62 defining the corrections of the transport mean free 
path in orders ln{k£) / {k£)^ and l/(fc£)^, see I I72I I. for the different contributions 
from the diagrams in figure [3] The data are obtained either analytically from I I64I I 
and I I67I I or by fitting from I I57I I or II65I I. respectively. 







SHA,2,t/i 




SBc+D,t/e 


5BE+F,t/i 














61 




1.234 


-2.968 


2.967 


3 


b2 


-(3/2 + 57r2/4) 


2.439 


0.906 


4.537 


2 - 51n(2) 



4-2. The transport mean free path 

As in the case of the scattering mean free path, the leading contribution originates 
from the diagram in figure U which is given by (j64p . The leading asymptotic behaviour 
of this contribution is SHA,i.t/i — > 57r/6fc^, for k£ — )■ 00. Like in the case of the 
scattering mean free path, this diagram also contributes to the logarithmic term 
(\n{k£) / {k£)^) and terms of order and higher. 

The remaining contributions are at least of logarithmic form and we thus write 
their asymptotic behaviour like in (|70l) : 

— —-^ {k£Y for 00. (72) 

Table [2] lists the coefficients 61 and 62 obtained either analytically (for (5i?E+F,t) or by 
fitting ([72]) on the numerically integrated contributions (cf. (|57|) and ([65])). 

Using (|55|) and collecting all terms up to order ln(fc£)/(fc€)^ we obtain the 
correction of the transport mean free path 

Similar to the case of the scattering mean free path (cf. end of section lCT , we analysed 

the behaviour of the weak localisation corrections due to long loops with at least three 

scattering events. This is achieved by inserting C2 = L^, instead of the full C, into 

(2) (2) (2) (2) 

the extended box H\ or the other boxes ^^b+c ' -^a+b '^"^^ ^c+d ' ^^'^ leads to 

vanishing logarithmic terms b\. As for the remaining coefficient 62, the contributions 
of the ladder boxes iJ^+B -^c+d cancel each other. In contrast, the Hikami boxes 

f21 (2) 

H\ 2 and ^^b+c '^o not compensate for each other and the resulting correction of the 
transport mean free path in this case reads Sltjl ~ —^Ajikt)^- This shows that the 
long loops contribute only to order l/(fc^)^, which is consistent with the standard 
result obtained within the diffusion approximation, see e.g. fP chapter 7.4 (7.57). 

A comparison of the Xjkl contribution in ((73)) with ([7T|) shows that 1^ is more 
strongly reduced than £f The relative increase of the transport mean free path can 
be made plausible by looking more closely at the leading diagram, see figure |31 and 
its effect on the step length distribution p6p . It turns out that the main contribution 
to this diagram comes from constellations where r\ and r2 are very close to each 
other. This approximately amounts to two independent steps, each of which follows 
the exponential step length distribution Pq (cf. (fT6|) ). at once. Due to this "double 
step" the variance of the step length distribution, and thus the transport mean free 
path, is increased with respect to the case with just single steps. 
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We finally look at the weak localisation corrections of the conductivity which were 
calculated in [13) : 

/ 27r TT^ -41n(H) 



(74) 

^0 I 1 - :=ri; - """la'"^' + 0((H)- ^ 



2. 0.7UHM) , 

ike {kef 



We compare this with the quantity V^^' — 4^t/3 for which we find, putting all 
previous corrections of i!s and together, the following leading order weak localisation 
corrections: 



-D''^ = ^ ( 1 " - ^"::zr' + omr') i . (75) 



1 27r 0.734 ln(H) ^ ^^_2^ 
3fc^ ^^)^ 



Hence, the leading weak disorder corrections of the conductivity a and of the stationary 
diffusion constant T)'^^^ are identical. Since, to the best of our knowledge, a general 
relation between I?*^^^ and a has not yet been proven, it is an interesting open question 
whether the proportionality between T)'^^^ and a remains valid also for stronger disorder 
(i.e. beyond the leading orders in Xjkt). 



5. Conclusion &: Outlook 



We calculated the leading weak disorder corrections of the scattering and of the 
transport mean free path for the propagation of waves in a white-noise random 
potential, without employing any approximation except fc^ 3> 1 (where e is the 
Boltzmann mean free path). The leading corrections for both, and i?t, scale like 
l/(fc£), and originate from the short-loop diagrams given in figure 2] In contrast, the 
self-consistent theory of localisation [181 El] , where the loops are treated within the 
diffusion approximation, predicts a vanishing correction of ^s, whereas the correction 
of 4 scales like Xjilitf . This is not necessarily a contradiction, since the self-consistent 
theory of localisation is meant to describe the transition from the regime of weak to 
strong localisation, rather than to calculate corrections for the case of weak disorder. 
The basic idea of the self-consistent theory is the renormalisation of the diffusion 
constant, where the loops contributing to the return probability (which, in turn, 
determines the diffusion constant) are treated in the diffusion approximation, but 
with a renormalized diffusion constant. Since the present paper shows that also short 
loops give rise to significant weak disorder corrections for both, £s and ^t, this may 
motivate future attempts to improve the self-consistent theory of localisation (which, 
in its present form, does not reproduce the correct critical exponents for the Anderson 
transition [231 124)). e.g. by renormalizing not only a single quantity describing 
transport on large length scales (such as the diffusion constant or the transport mean 
free path), but taking into account also quantities relevant for transport on shorter 
length scales (such as the scattering mean free path). 

Furthermore, let us comment on the relevance of our results with respect to 
experiments on weak localisation [7ll25[. The experiments described in T consider the 
weak localisation corrections of the electron conductivity (or, similarly, the mobility) 
under breaking of the reciprocity symmetry, whereas the disorder parameter l/fc£ is 
kept constant. Breaking of the reciprocity symmetry is achieved by the application of 
a magnetic field which affects mainly the long interference loops, since these encircle 
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on average large areas being sensitive to penetration by the magnetic flux. Thus, these 
experiments observe the weak localisation contribution of order characteristic 
for long loops (the logarithmic correction being hard to distinguish from the pure 
l/{ki)'^ decay, experimentally). However, we showed that the diagram depicted in 
figure [4] gives rise to a correction of the order 1 / k£. Again, this is no contradiction 
since this diagram is short range and thus almost insensitive to reciprocity breaking 
induced by the magnetic field. Indeed, other experiments [26], which measure the 
mobility of electrons in helium gas and directly change the density of the scattering 
medium and thus the disorder parameter 1/ki, were confirmed to be consistent with 
a leading order correction to scale like l/k£ in [25] . 

Based on the thorough understanding of the linear provided in the 

present article, it becomes now furthermore possible to investigate the influence of non- 
linearities on the weak localisation effect. Application of the non-linear diagrammatic 
theory developed in [^71 [55] to the case of non-linear coherent backscattering showed 
an inversion of the CBS cone, already for moderate values of the non-linearity [291 . 
This inversion stems from destructive interference in the backscattering direction and 
thus suggests a weak anti- localisation effect. This cone inversion was also observed in 
numerical experiments on coherent backscattering of Bose- Einstein condensates in two 
dimensional speckle potentials |30| . in accordance with the non- linear diagrammatic 
theory. Combining the non-linear diagrammatic theory and the work presented in this 
article it is possible to calculate the corrections of the transport mean free path in the 
non-linear case. The negative cone should then correspond to an increased transport 
mean free path as compared to the classical (ladder approximation) prediction. 
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